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ABSTRACT 

The infall regions of galaxy clusters represent the largest gravitationally bound structures in a 
ACDM universe. Measuring cluster mass profiles into the infall regions provides an estimate of the ul- 
timate mass of these haloes. We use the caustic technique to measure cluster mass profiles from galaxy 
redshifts obtained with the Hectospec Cluster Survey (HeCS), an extensive spectroscopic survey of 
galaxy clusters with MMT/Hectospec. We survey 58 clusters selected by X-ray flux at 0.1<z<0.3. 
The survey includes 21,314 unique MMT/Hectospec redshifts for individual galaxies; 10,275 of these 
galaxies are cluster members. For each cluster we acquired high signal-to-noise spectra for ~ 200 
cluster members and a comparable number of foreground/background galaxies. The cluster members 
trace out infall patterns around the clusters. The members define a very narrow red sequence. The 
velocity dispersions decline with radius; we demonstrate that the determination of the velocity disper- 
sion is insensitive to the inclusion of bluer members (a small fraction of the cluster population) . We 
apply the caustic technique to define membership and estimate the mass profiles to large radii. The 
ultimate halo mass of clusters (the mass that remains bound in the far future of a ACDM universe) 
is on average (1.99±0.11)M2oo, a new observational cosmological test in essential agreement with 
simulations. Summed profiles binned in M 2 oo and in Lx demonstrate that the predicted NFW form 
of the density profile is a remarkably good representation of the data in agreement with weak lensing 
results extending to large radius. The concentration of these summed profiles is also consistent with 
theoretical predictions. 

Subject headings: galaxies: clusters: individual — galaxies: kinematics and dynamics — cosmology: 
observations 



1. INTRODUCTION 

Clusters of galaxies are the most massive virialized sys- 
tems in the universe. Clusters are surrounded by infall 
regions where galaxies are bound to the cluster but are 
not in dynamical equilibrium in the cluster potential. If 
dark energy behaves like a cosmological constant, clus- 
ter infall regions are the largest gravitationally bound 
structures in the universe. Thus, measuring cluster mass 
profiles at large radii provides an estimate of the ulti 



mate halo mass of these systems jNagamine fc Loeb 2003 
Busha et al.||2005| |Diinner et al [|2U0B ) 
The Cluster and" Infall Region 



Nearby Survey 
(CAIRNS) pioneered the detailed observational study 
of cluster infall regions. CAIRNS studied nine nearby 
galaxy cluste rs and their infal l regio ns with extensive 
spectroscopy ( Rines et aT]|2003| |2005[ ) and near-in frared 
photometry from the Two-Micron All-Sky Survey (Rines 
et aH2004j ). The nine CAIRNS clusters display a charac- 
teristic trumpet-shaped pattern in radius-redshift phase 
space diagrams. These patterns wer e first predicted for 



simple spherica l infall onto clusters (Kaiser 

1 



1987 Regos 



fe Geller] 1989 ) , but later work showed that t hese pat 
terns reflect the dynamics of the infall region (Diaferio 
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fc Geller|p97] |Diaferio|p99l hereafter DG and D99). 
The Cluster In fall Regi ons m the Sloan Digital Sky Sur- 
vey ( |Rines fc Diaferio||2006| hereafter CIRS) project ex- 
tended this analysis to 72 X-ray -selected clusters in th e 
Sloan Digital Sky Survey (SDSS, |Stoughton et al |2002 ). 
CIRS showed that these infall patterns are ubiquitous in 
nearby X-ray clusters. 

Using numerical simulations, DG and D99 showed that 
the amplitude of the caustics is a measure of the escape 
velocity from the cluster; identification of the caustics 
therefore allows a determination of the mass profile of 
the cluster on scales < 10/i _1 Mpc. In particular, non- 
parametric measurements of caustics yield cluster mass 
profiles accurate to ~50% on scales < 10/i — 1 Mpc when 
applied to Com a-size clusters extr acted from cosmologi- 
cal simulations. ISerra et al. (2011) confirm these results 



for clusters across a broader mass range and they show 
that the dominant source of uncertainty in individual 
cluster mass profiles is projection effects from departures 
from spherical symmetry. The caustic technique assumes 
only that galaxies trace the velocity field. Indeed, Simula- 



and mildly non- 


linear scales (Kauffmann et al.| 


1999a b 


Diemand et al. 


2004 


Faltenbacher et al. 200c 


1 This 



conclusion is supported observationally by the excellent 
agreement between the clu ster virial mass fu nction and 
other cosmological probes ( Rines et al.||2008 ). 

CAIRNS and CIRS showed that caustic masses of clus- 
ters agree well with mass estimates from bot h X-ray ob- 
serva tions an d Jeans' analysis a t smal l radii (Rines et al. 
20031 CIRS). |Lokas fc Mamon| (|2003| confirm that the 
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Rines et al. 



mass of Coma estimated from higher moments of the 
velocity distribution agrees well with the caustic mass 
estimate ( |Geller et aL|l999 ). 

The caustic technique provides an estimate of the mass 
profile of clusters. For instance, CAIRNS and CIRS 
showed that cluster mass profiles fall off more steeply at 
large radii than an isothermal sphere. Thus, caustic mass 
profiles probe the structure of dark matter haloes, and 
these profiles can be comp ared to those determ ined from 
gravitational lensing (e.g., Umetsu et al.|[20lT Geller et 
al. 2012). 

At large radii, neither galaxies nor intracluster gas 
should be in equilibrium, thus invalidating the use of the 
virial theorem or hydrostatic equilibrium at these radii. 
The caustic technique and gravitational lensing are the 
only cluster mass estimators that do not rely on the equi- 
librium assumption. Gravitational lensing is necessarily 
contaminated by linc-of-sight structure unrelated to the 
clust er; this contamination becomes larger at larger radii 
(e.g., [Hockstra et al. 201^0 . Despite this potential diffi- 
culty, Dia ferio et al.| ( |2005p showed that caustic masses 
agree with weak lensing masses in three clusters at mod- 
erate redshift. 

CAIRNS and CIRS showed that infall patterns are well 
defined in observations of nearby massive clusters. In 
fact, the infall patterns or "caustics" have significantly 
higher contrast in the CAIRNS observations than in the 
simulations of DG and D99. The CAIRNS and CIRS 
clusters are fairly massive clusters and gener ally have lit- 



tle surroun ding large-scale structure (but see Rines et al. 



2001 2002 



One might suspect that the presence of infall patterns 
is limited to massive, isolated clusters. However, other 
investigators have found infall pattern s around the For 



nax Cluster (Drinkwater et al. 



percluster (fReiscncgger et al. 



2001) 



the Shapley Su- 
an ensemble clus- 
ter comprised of poor cl usters in the Two Degree Field 



Galaxy Redshift Survey (Biviano & Girardi 2003 ), and 



even the galaxy group associated with JN GC 5846 ( Mah 
ivi et al. 2005[ ). Caustics are also iden tifiable in X-ray 
groups m SDSS dRines & Diaferio||2010|). 

The presence ot caustics m all nine CAIRNS clusters 
and in all 72 CIRS clusters suggests that they are ubiq- 
uitous in nearby, massive clusters. Because clusters and 
especially cluster infall regions form quite late in the evo- 
lution of structure, infall patter ns evolve even at m o dest 
redshifts (se e an a nimation in Geller et al. 20111. Di- 
aferio et al.| ( |2005[ ) showed that intall patterns exist in 
three clusters at moderate redshift, althou gh these clus- 
ters w ere not carefully selected. Similarly, Lemze et al. 
(2009) showed that the caustic mass profile of A1689 



(z=0.18) agrees with the mass profile determined from a 
joint analysis of X-ray and lensing data. 

Inspired by these results, we conducted a systematic 
survey of infall regions at z«0.2. The Hectospec Clus- 
ter Survey (HeCS), includes MMT/Hectospec spectra of 
large samples of galaxies in the infall regions of X-ray- 
selected clusters at z=0. 1-0.3. HeCS is the first sys- 
tematic spectroscopic survey of cluster infall regions at 
z>0.1. 

Spectroscopic data for clusters at moderate redshift 
can test the robustness of scaling relations between differ- 
ent cluster observables that correlate with cluster mass. 
For instance, the integrated Sunyaev-Zeldovich decre- 



ment Ysz has small scatter in simulations (Hallman 
et al.|200"6 ), as does th e X-ray observable Yx =~M, 



J-x 



\ Kravtso'v'et al.||2006 1 . The Yx parameter also seems to 



have low scatter in observations, although the esti mated 
masses and Yx are derived fro m the same data (Nagai 
|et al.||2007| ). |Rines et ah ( 2010[ ) studied the relation be- 
tween Ysz and the dynamics of galaxies in a sample of 
15 clusters. The dynamical properties correlate strongly 
with increasing Ysz but with significant scatter. 

We describe the data and the cluster sample in § 2. 
In § 3, we review the caustic technique, use it to es- 
timate the cluster mass profiles, and estimate the ulti- 
mate halo masses of clusters. To mitigate systematic 
uncertainties due to projection effects, we analyze en- 
semble clusters in § 4 and compare the ensemble mass 
profil es with the theoretical NFW profile (|Navarro et al. 



19971 and with profiles determined from gravitational 
lensing. We discuss the color distribution of cluster 
members in § 5 and show that our red-sequence selec- 
tion does not bias our dynamical measurements. We 
discuss our results and conclude in §6. We assume 
H = 100ft km s _1 ,fi m = 0.3, fl A = 0.7 throughout. We 
measure cluster masses Ma , defined as the mass enclosed 
within the radius r/\ that encloses an average density of 
Ap c (z), where p c {z) = 3H 2 (z) /8irG is the critical den- 
sity at redshift z. 

2. THE HECS CLUSTER SAMPLE 

We construct the HeCS sample to take advantage of 
two large-a rea public surveys: the Sloan Digital Sky Sur- 
vey (SDSS, |Stoughton_efal][20^1 iand the ROSAT All- 



Sky Survey (RASS, |Voges et al.||1999| . In particular, 
we utilize existing X-ray cluster catalogs based on RASS 
data to define a flux-limited cluster sample. We then 
match these systems to the imaging footprint of the SPS S 
Data Release 6 (DR6, |Adelman-McCarthy et al.||2008| . 
The accurate SDSS multicolor photometry enables selec- 
tion of candidat e cluster members usin g the red-sequence 
technique (e.g., Gladders & Yee 20001. 

We obtained MlVLT/ Hectospec spectroscopy of 400-550 
candidate members per cluster. Hectospec redshifts en- 
able robust membership classification and estimates of 
the virial masses of the clusters. The wide field-of-view 
of Hectospec allows us to simultaneously probe the virial 
and infall regions of these clusters. 

The earlier CIRS project used spectr oscopic data from 
the S loan Digital Sky Survey (SPSS IStoughton et a~ 
20021 to study clusters at z<0.1 identified in X-ray clus- 
ter catalogs based on the ROSAT All-Sky Survey. At 
z>0.1, the SPSS redshift survey is not dense enough for 
accurate assessment of cluster masses. 

2.1. Sloan Digital Sky Survey 



Th e Sloan Pigital Sky Survey (SPSS, |Stoughton et al. 
2002 ) is a wide-area photometric and spectroscopic sur- 
vey at high Galactic latitudes. The Sixth Pata Release 
(PR6) of SPSS includes 8417 square degrees of imaging 
data ( j Adelman- McCarthy et al ||2008[ ) . 

From a comparison of SPSS with the Millennium 



Galaxy Catalogue, Cross et al. ( 2004 ) conclude that there 
is a photometric incompleteness ot ~7% due to galaxies 
misclassified as stars or otherwise missed by the SPSS 
photometric pipeline. For our purposes, the incomplete- 



HcCS: Ultimate Halo Mass 



3 



ness is not important provided sufficient numbers of clus- 
ter galaxies do have spectra. Further, the cluster galax- 
ies we focus on here are intrinsically luminous and have 
higher surface brightnesses than less luminous galaxies. 
Thus, the photometric incompleteness is probably some- 
what smaller for these galaxies. 

The spectroscopic limit of the main galaxy sample of 
SPSS is r=17.77 aft er correcting for Galactic extinction 
([Strauss et al.||2002| . CAIRNS and CIRS show that in- 
tall patterns are easily detectable in clusters sampled to 
about M* + 1, or z <0.1 for the SDSS spectroscopic 
limit. Studying cluster infall regions at z~0.2 requires 
spectroscopy to a limit of r ~ 21. 

2.2. X-ray Cluster Surveys 
We select the HeCS clusters from X- ray catalogs based 



2.3. MMT/Hectospec Spectroscopy 



on the ROSAT All-Sky Survey (RASS, |Voges et al. 1999| ). 
RASS is a shallow survey, but it is sufficiently deep to 
include nearby, massive clusters. RASS covers virtually 
the entire sky and is thus the most complete X-ray cluster 
survey for nearby clusters. The flux limits of RASS-based 
surveys are fti3xl0~ 12 erg s -1 cm~ 2 in the RO SAT band 
( |Ebeling et al.||1998| [Bohringer et aT|2000b|a| . 

We restrict our study to systems with 0.10<;z<0.30. 
The flux-limited HcCS sample consists of the 53 clusters 
from the BCS and REFLEX catalogs within the SDSS 
DR6 photometric footprint and with /x>5x 10 _12 erg 
s _1 cm -2 . Due to scheduling constraints on observing 
time, HeCS includes four additional clusters (A2187, 
A2396, A2631, and A2645) with / x >3x 10~ 12 erg 
s cm 



A2645 is in REF LEX and the others a re from 
the extended BCS (eBCS; [Ebeling et aT1|2"000[). BCS 



splits A1758 into two components, neither of which alone 
would lie above our flux limit ; however, the NORAS cat- 
alog ( |Bohringer et al. 2000b[ ) merges these components, 
yielding a ffux above our flux limit. We also include 
an X-ray cluster, A750 , that lies in the foreground of 



MS0906+H10 (see |3_6)for details). HcCS thus contains 
58 clusters: 53 clusters in the flux-limited sample, four 
clusters with slightly smaller X-ray fluxes, and A750 (we 
do not count the subclusters of A1758 as individual clus- 
ters). 

A689 has a larg e contribution from a central BL Lac; 



Giles et al.| ( |2012[ ) estimate that only 10% of the BCS 
luminosity is from the intracluster medium. The long 
arrow in Figure [T] shows that this cluster would not lie 
in the flux-limited sample with the corrected luminos- 
ity. Similarly, the luminosity of A1758N (the brighter of 
A1758N/S) alone would remove it from the flux-limited 
sample (short arrow in Figure [lj. 

Table [l] describes the basic properties of the HeCS clus- 
ters. Figure n] shows the (rest-frame) X-ray luminosities 
of the HeCS clusters compar ed to the HiFluGCS survey 
( Reiprich fc B6hringer|[2002| and to our previous survey 
ClRS. HeCS contains systems with significantly larger 
Lx than CIRS, although there is a substantial over- 
lap at intermediate Lx, especially at 2=0.10-0.15. The 
larger redshift limit of HeCS relative to CIRS probes a 
much larger volume: the HeCS volume is ^10 8 /i~ 3 Mpc 3 , 
a factor of >10 larger than the volume probed by CIRS. 
The number density of clusters declines with Lx (e.g., 
|Bohringer et al. 2002); a larger survey volume increases 
the sample of intrinsically more luminous systems. 



The Hectospec instrument ( Fabricant et al. 2005 1 
mounted on the MMT 6.5m telescope is ideal for study- 
ing cluster infall regions at moderate redshift. Hectospec 
is a multiobject fiber-fed spectrograph with 300 fibers 
deployable over a circular field-of-view with a diameter 
of 1°. One Hectospec pointing extends to a radius of 
2.3 (5.6) h~ 1 Mpc at z=0.1 (0.3), so a single Hectospec 
pointing covers the entire virial region and extends well 
into the infall region. 

Our observing strategy for HeCS was to identify the 
red sequence of cluster galaxies for each system and tar- 
get primarily galaxies withi n 0.3 magnitudes of the red 
sequence. As shown in §5.1[ the color selection is signifi- 
cantly broader than the actual red sequence of confirmed 
cluster members. Thus, the target selection includes pri- 
marily, but not exclusively, red-sequence galaxies. We 
use SDSS DR6 spectroscopic data to identify galaxies 
with existing redshifts and remove them from the target 
catalog. The galaxy targets have r magnitudes of r=16- 
21. The Hectospec fiber assignment software xfitfiba^ al- 
lows the user to rank targets with priorities. We ranked 
galaxies primarily by their proximity to the cluster center 
and secondarily by apparent magnitude. This procedure 
yields a largely complete magnitude-limited sample of 
brighter galaxies supplemented by a more sparsely sam- 
pled selection of galaxies up to 1 mag fainter than the 
bright limit. The Appendix describes the target selec- 
tion procedure in more detail. 

The galaxy targets are relatively bright for Hectospec 
spectroscopy: we can obtain high-quality spectra with 
3x20-minute exposures even under suboptimal observing 
conditions (e.g., poor seeing, thin clouds). Because Hec- 
tospec is a queue-scheduled instrument, this flexibility 
allows the HcCS fields to be observed under many ob- 
serving conditions and improve the overall efficiency of 
the Hectospec queue. Observations were conducted pri- 
marily between 2007 June and 2009 February with a total 
of 10 nights of queue time. We also observed one addi- 
tional Hectospec field in A2261 in 2011 May and seven 
additional fields in RXJ2129 in 2011 September. We ob- 
served additional fields in the se clusters because they are 
part of the CLASH sample (Postman et al. 2012 Coe 



et al.||2012[ ). The supplementary Hectospec fields tar 
geted galaxies of all colors, enabling us to test the poten- 
tial impact of our red-seque nce selection on the estimates 



of dynamical parameters Q5.3). 

Because the number of galaxies per cluster varies sig- 
nificantly, we adjusted the limiting magnitudes primarily 
to obtain large, nearly complete, samples to the faintest 
magnitude possible. Thus, the limiting absolute magni- 
tude varies significantly from cluster to cluster. 
After processing and r educing the spectra, we used 



1998) to cross- 



the IRAF package rvsao (Kurtz & Mink 
correlate the spectra witn galaxy templates assembled 
from previous Hectospec observations. During the 
pipeline processing, spectral fits are assigned a quality 
flag of "Q" for high-quality redshifts, "?" for marginal 
cases, and "X" for poor fits. Repeat observations of sev- 
eral targets with "?" flags show that these redshifts are 
generally reliable. We visually inspected all spectra to 



5 https : //www. cf a. harvard. edu/~john/xf itf ibs/ 
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HeCS Basic Properties 



Name 


X-ray Coordinates 




z® 




Catalog 


CTp 






ra n2oooi 


DEC ('12000') 






erg s ^ 




km s ^ 




A 9R7 


28.1762 


01.0125 


n 
u 


99Q1 

.zzy i 


A 1 
10 




y '^-53 


99R 
ZZO 


"7t,t1 47Q 


119.9190 


53.9990 


n 

u. 


1 097 


O flA 


t-> O o 


47Q+ 66 
4 ' y -46 


S9 
OZ 


A R4R 


125.5470 


47.1000 


u 


1 973 
.1Z I O 


1 99 
l.ZZ 




boci_ 51 


9(^4 
ZO^t 


A R 1 ^ 
/\000 


126.3610 


47.1320 


n 

u. 


1 971 


1 Q1 




777+58 
' ' '-47 


31 ^ 
OlO 


A RR7 


127.0190 


44.7640 


u 


1 4^9 

. l^OZ 


1 39 
1 .oz 




b4o_ 58 


1 4& 


A R£Q 


129.3560 


14.9830 


n 

u. 


97SQ 

.z / oy 


Q f!9 

y .oz 




o»y_ 62 


1 R3 
IOO 


A RQ7 


130.7362 


36.3625 


n 

KJ- 


981 9 
.ZolZ 


^ 1 ^ 




1 009~l~ 9 ^ 
lUUZ_yc [ 


iok 
loO 


A 7^0 


137.2469 


11.0444 


n 
u. 


1 R40 
. lO^tU 




Dpq* 


0»1_ 45 


99^ 
ZZO 


lVfCIOQOR 

iviouyuo 


137.2832 


10.9925 


n 
KJ 


1 7R7 


3 93 
O.ZO 




bb4_ 62 


1 01 
1U1 


A 773 


139.4624 


51.7248 


n 
u. 


91 73 
Zl ( O 


3 OS 

o.yo 


±3 O O 


iiiu_ 70 


1 73 
1 1 


A 70^ 


141.0240 


14.1680 


n 
u 


1 374 


1 

1 .oo 




7VQ+61 

' '^-50 


1 70 

1 ( y 


7w9701 


148.1980 


51.8910 


u. 


91 RO 
ZIOU 


3 30 
O.OU 




OOZ 55 


03 

yo 


A QR3 


154.2600 


39.0484 


u 


9041 


3 07 
O.U i 




0^+ 80 
yob_ 64 


91 1 

Zll 


a oso 


155.6275 


50.1017 


KJ. 




9 OR 
Z.UO 




1 033"'"'''^ 

lUOO_r,g 


999 
zzz 


ZjWOI^O 


155.9117 


04.1865 


n 

u. 


98Q4 


& 38 
o.oo 




op-o+103 


1 OR 
1UO 


A QQO 

Ayyu 


155.9120 


49.1450 


n 
u. 


1 41 R 

. 1^10 


9 13 
Z. IO 


DpC 


boo_ 60 


01 

y 1 


7 W 31 7Q 


156.4840 


12.6910 


n 
KJ 


1 499 


1 34 
1 .O^t 




r 41 +122 
041_ 73 


oy 


A i 033 


157.9320 


35.0580 


n 
u. 


1 990 
. 1ZZU 


1 3^ 
1 .OO 


13 O O 


^77+55 
b ( / _44 


1 01 

lyi 


A i or& 


160.1870 


39.9510 


n 
U 


1 3£R 


O 1Q 
Z. lo 


T3P^ 


109S + 106 


1 90 

izy 


A 1 1 39 


164.6160 


56.7820 


u. 




1 09 

i .yz 


13 O O 


' 4y -61 


1 RO 
lOU 


A 1 901 


168.2287 


13.4448 


u 


1 R71 
. ID ( 1 


1 70 

i . i y 


RPd 
13 V_ O 


fi oo+68 


100 


A 1 904 


168.3324 


17.5937 


KJ. 


1 70R 
. 1 / UO 


9 13 
Z. IO 


Rpq 

13 »w O 


^qo+62 
ooz_ 4g 


09 

yz 


A 1 93^ 


170.8040 


19.6160 


n 
u 


1 030 


O R^ 

u.oo 


T3P^ 

13 O O 


co 4 +62 
004_ 47 


1 31 
101 


A 1 94R 

.TilZ'iO 


170.9912 


21.4903 


n 

u. 


1 Q91 

. iyzi 


9 31 
Z.Ol 


Rpq 

13 O O 


Q0R+ 70 

yuo_5 7 


99R 
ZZO 


A 1 309 


173.3070 


66.3990 


n 
U 


1 1 ^9 
. 1 lOZ 


O fiA 


RfC 
13 O O 


D0U _48 


1 R9 
lOZ 


A 1 3R1 


175.9170 


46.3740 


n 
u. 


1 1 

. i ioy 


O QQ 

u.yy 


Rpq 

13 »w O 


k-io+64 
Ol^_ 47 


iyo 


A 1 3RR 


176.2020 


67.4130 


u 


1 1 £0 


1 ns 

1 .Uo 


Rpq 

13 O O 


fi1 «+62 


900 

zuu 


A 1 41 3 


178.8260 


23.4080 


n 

u. 


1419 


3 71 
O. ( 1 


Rpq 

13 O O 


OCfi+90 

»ob_ 68 


1 1 fi 
110 


A 1 493 


179.3420 


33.6320 


u 


91 49 


3 07 
O.U / 


RPC 


7C rQ+64 


930 

zou 


A 1 437 


180.1040 


03.3490 


KJ- 


1 333 

. lOOO 


9 19 
Z. 1Z 


Rpq 

13 »w O 


IZOJg^ 


1 04 

ly^i 


A 1 ^3 


187.6959 


10.5606 


n 
u. 


■ lOOO 


9 17 
Z. 1 ( 


Rpq 

13 O O 


8b '-51 


1 71 
1 1 1 


A 1 ££9 


196.7278 


46.5560 


u. 


9979 
,ZZ 1 z 


3 4S 


RPC 
13 O O 


yyb -65 


1 ^1 
101 


A 1 R8Q 


197.8750 


-1.3353 


n 
U 


1 S49 


7 OR 
/ .UO 


ixUjI' 


119 '_65 


910 
Z1U 


A 1 7^8 


203.1796 


50.5496 


n 
KJ. 


97(^0 
Z / OU 


^ 89 
O.oZ 


RfC 
13 O O 


R74+" 


1 43 
I'iO 


A 1 7R3 


203.8257 


40.9970 


u. 


931 9 
.ZolZ 


4 79 


T3P^ 


i np. i +81 
1261_ 6g 


937 
ZO / 


A 1 £3^ 


210.2595 


02.8801 


n 

u. 


ZOUO 


1 1 77 
11. ( ( 


RPC 
13 »w O 


iicri+80 
llt)l_ 66 


910 

ziy 


A 1 009 


215.4226 


37.2958 


u. 


1 fi93 
. lOZO 


1 .0 i 


RPQ 
13 O O 


' 84 -56 


1 30 
lOU 


A 1 Q1 8 


216.3420 


63.1830 


n 
u. 


1 3£8 


117 
1.1 ( 


Rpq 

13 O O 


04t)_ 54 


so 

0U 


A 1 Q1 4 


216.5068 


37.8271 


n 
u. 


. lOOU 


^ 03 
O.UO 


RfC 




9^ 
ZOO 


A i G30 


218.1200 


31.6330 


n 
u. 


1 30& 


11^ 
1.10 


Rpq 

13 O O 


c 77 +75 
0' '-54 


7fi 
1 


A 1 Q7& 

/\iy f o 


222.7750 


14.6110 


n 
U 


1 4^Q 


1 90 

i .zy 


Rpq 

13^ O 


4U4 -56 


R3 
OO 


A 9000 


225.0850 


21.3620 


n 
u. 


1 ^99 
. 10ZZ 


9 fiO 
Z.OU 


Rpq 

13 »w O 


7-1C+57 


iyo 


T?"V T1 ^04 


226.0321 


-2.8050 


u. 


91 I^S 


14 19 
1*±. 1Z 


DTP 

ixUjI' 


77O+105 
' ' y -75 


1 90 

1ZU 


A 9034 


227.5450 


33.5060 


n 

u. 


1 1 39 
. 1 lOZ 


1 09 

i .yz 


Rpq 

13 »w O 


q 49 +64 
y4 ^-53 


1 89 
loZ 


A 90^0 


229.0680 


00.0890 


n 
u. 


1 1 Q1 

. i iy i 


1 IO 

1.1(5 


RPC 
13^ O 


Qfio+77 


1 OR 
1U0 


A 90^ 


229.6720 


06.2110 


u. 


1 093 
. 1UZO 


1 34 


Rpq 

13 O O 


fi7fi+ 90 
b ' b -64 


930 
ZOU 


A 90RQ 


231.0410 


29.9210 


KJ 


1 1 3Q 

. i ioy 


9 4fi 


Rpq 

13 »w O 


yy4 -52 


441 


A 91 1 1 
iAZl 1 1 


234.9337 


34.4156 


n 

u. 


99Q1 

zzy i 


3 3^ 
O.OO 


RPC 
13 O O 


7A1 + 65 
' 4i -52 


90S 
ZUo 


A 91 £7 


246.0591 


41.2383 


n 
U 


1 &9Q 

. lozy 


1 ^fi 
1.00 




rqi +83 


1 03 
1UO 


A2219 


250.0892 


46.7058 


o 


.2257 


6.10 


BCS 


1151+S 


461 


Zw8197 


259.5480 


56.6710 





.1132 


0.80 


BCS 




76 


A2259 


260.0370 


27.6702 


0. 


.1605 


1.85 


BCS 


8551™ 


165 


RXJ1720 


260.0370 


26.6350 





.1604 


4.47 


BCS 


860t« 


376 


A2261 


260.6129 


32.1338 





,2242 


5.55 


BCS 


780t^ 


209 


RXJ2129 


322.4186 


00.0973 





.2339 


5.65 


BCS 


858+^ 


325 


A2396 


328.9198 


12.5336 


0. 


.1919 


1.86 


eBCS 


935+?;; 


176 


A2631 


354.4206 


00.2760 





.2765 


4.15 


eBCS 


851+ 7 « 


173 


A2645 


355.3200 


-9.0275 





.2509 


2.85 


REF 


549+^ 


61 



HcCS: Ultimate Halo Mass 



■5 



include "?" and "X" spectra that have redshifts secured 
by multiple lines (usually four or more). Repeat observa- 
tions indicate that the redshift uncertainties from rvsao 



are reasonable: |Geller et al. (2012 1 use 1468 unique pairs 
of repeat observations to estimate a mean internal error 
of 56km s _1 for absorption-lin e objects and 21km s - 1 for 



emission-line objects (see also Fabricant et al. [20051 
Here we include new redshifts for 57 of the HeC 



ters. For th e remaining cluster, RXJ1720- 



eC S clus- 
26, |Owers| 

et al I ( |2011[ > obtained redshifts from MMT/Hectospec 
and Keck/DEIMOS for a separate investigation of cool- 
core clusters. The data for RXJ1720+26 do not extend 
to the full field of Hectospec, but the spectra do extend 
to fainter apparent magnitudes than HeCS. 

We observed clusters at z>0.15 with two configurations 
to obtain larger samples and to mitigate issues with fiber 
collisions. We observed clusters at 2=0.10-0.15 with one 
Hectospec configuration. Table [2] lists the 21,314 unique 
redshifts obtained for HeCS. In our analysis, we also use 
~7,000 reds hifts from SPSS an d 132 redshifts around 
A2219 fromlBoschin et al.l (|2004h. 



3. RESULTS 
3.1. Ubiquity of Infall Patterns around Clusters 

We first search for well-defined infall patterns around 
X-ray clusters. Analogous to CIRS, we plot the rest 
frame line-of-sight velocity relative to the cluster cen- 
ter as a function of projected radius in Figure [2j All 58 
HeCS systems have "clean" infall patterns; that is, there 
is little ambiguity in the location of the caustics or limits 
of the pattern in redshift. All clusters contain a large 
number of members at the cluster redshift extending out 
to several Mpc from the cluster center. 

Figure[T]shows the X-ray l uminosities of the HeCS clus- 
ters compared to HIFluGCS ( |Reiprich fc Bohringer|2"002 1 
and CIRS. Due to the larger survey volume of HeCS, it 
includes many more systems with Lx^l0 44 /i~ 2 crg s _1 
than the CIRS sample. There is a substantial overlap at 
intermediate Lx, especially at z=0.10-0.15. Table [T] lists 
the clusters in the HeCS sample, their X-ray positions 
and luminosities, their central redshifts and rest-frame 
velocity dispersions. Velocity dispersions are measured 
for galaxies within r2oo as determined from the caustic 
mass profiles. The radius ta is the radius within which 
the enclosed average mass density is Ap c , where p c is the 
critical density) by computing the enclosed density pro- 
file [p(< r) = 3M(< r)/4irr 3 ]; r2oo is the radius which 
satisfies p(< r 2 oo) = 200p c . 

Figure [2] shows the infall patterns and caustics for 
the HeCSclusters. The contrast in phase space density 
between cluster members and foreground/background 
galaxies is striking. 

The clusters are ordered by decreasing X-ray luminos- 
ity. Velocity dispersions of clusters can be inferred visu- 
ally from the spread in velocities at small radius. One 
immediate conclusion is that velocity dispersion and X- 
ray luminosity are not perfectly correlated: clusters such 
as RXJ1504 and A689 have large X-ray luminosities, but 
their velocity dispersions are smaller than many clusters 
with comparable X-ray luminosity. 

3.2. Caustics and Mass Profiles 




Fig. 1. — Redshift versus X-ray luminosity (0.1-2.4 keV) for X- 
ray clusters from CIRS (small open squares) and HeCS clusters 
(filled blue circles) contained in the SDSS DR6 imaging survey 
region. The blue solid lines show the selection of the HeCS sample: 
redshifts 0.1<z<0. 3 and a flux limit of f x >5x 10~ 12 erg s _1 . Four 
HeCS clusters are from a subsample of clusters with HA>17 h and 
flux .fx >3xl0~ 12 ergs- 1 (the same flux limit as CIRS). The lower 
solid line shows the flux and redshift limits of the CIRS cluster 
sample. The dash-dotted line shows the flux limit (2 xl0 _11 crg 
s" 1 ) of the HiFluGCS sample. 



We calculate the shapes of the caustics with the tech- 
nique described in D99 using a smoothing parameter of 
q = 25. The smoothing parameter q is the scaling be- 
tween the velocity smoothing and the radial smoothing 
in the adaptive kernel estimate of the underlying phase 
space distribution. Previous investigations show that the 
mass profiles are insensiti ve to changes of a factor of 2 in 
the smooth ing parameter ( Geller et al.|1999 Rines et al 
2000, p002 ). 

The technique of D99 uses the redshifts and coordi- 
nates of the galaxies to determine a hierarchical center 
based on a binary tree analysis. Analysis of 3000 simu- 
lated clusters indicates that the binary tree analysis re- 
covers the input cluster center to within ■^300/z ~ 1 kpc 
for 95% of simulated clusters (Serra et al. 2011). For 
the HeCS clusters, the hierarchical centers are located 
within 300/i _1 kpc of the X-ray coordinates for all but 
four clusters (7% of the sample, consistent with the 5% 
of simulated clusters with offsets this large). We discuss 
these four clusters and some other individual cases in 
For one cluster (A2261), we use the slightly dif- 



^ 3.6 



ferent algor i thm fo r cutting the binary tree described by 
Serra et al. (2011) to determine the center. 

Figure [2] shows the caustics. The D99 algorithm we use 
to identify the caustics generally agrees with the lines one 
would draw based on a visual impression. This consis- 
tency suggests that systematic uncertainties in the caus- 
tic technique are dominated by proje ction effects rathe r 
than the details of the algorithm (see Serra et al. 2011 \. 

Figure [3] shows the associated caustic mass profiles. In 
redshift space, a cluster of galaxies app ears as trumpet- 
shaped pattern ( |Regos fc Geller||1989[ ). DC and D99 
demonstrated that tor clusters forming hierarchically, the 
boundaries of this sharply defined pattern (termed caus- 
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TABLE 2 
HeCS Redshifts 



Coordinates 


CZQ 




R 


Flag 


RA (J2000) 


DEC (J2000) 


km s 


km s _1 






1:51:35.911 


0:37:43.140 


114091.06 


23.75 


16.22 


Q 


1:51:01.337 


0:51:26.568 


107171.56 


39.25 


7.61 


Q 


1:52:35.819 


0:55:05.660 


18329.14 


16.37 


18.88 


Q 


1:51:30.070 


0:52:59.772 


77724.84 


28.48 


18.52 


Q 


1:52:22.812 


0:52:20.964 


69890.66 


48.43 


8.65 


Q 



a Table[2]is published in its entirety in the electronic edition of the Journal. A portion is shown here for guidance regarding its form and 
content. 



tics) in redshift space (phase space) can be identified with 
the escape velocity from the cluster. This identification 
provides a route to estimation of the cluster mass pro- 
file assuming spherical symmetry. This mass estimation 
method is called the caustic method. 

The amplitude of the caustics A(r) is half the distance 
between the boundaries of the cluster in redshift space. 
With the assumption of spherical symmetry the gravita- 
tional potential <p(r) and the caustic amplitude A(r) are 
related by 



A 2 (r) = -20(r) 



1 - g(r) 
3 - 2p(r) 



(1) 



where /3(r) is the anisotropy parameter, /3(r) = 1 — 
Og(r)/[2af(r)] where ag and ay are, respectively, the tan- 
gential and radial velocity dispersions. 

DG show that the mass of a spherical shell within the 
infall region is the integral of the square of the caustic 
amplitude A(r): 



GM(<r)-GM(<r ) =T P f A 2 (x)da 

J rn 



(2) 



where J-p ~ 0.5 is a filling factor with a value estimated 
from numerical simulations. We approximate Tp as a 
constant; variations in J-p with radius lead to some sys- 
tematic uncertainty in the mass profile we derive from the 
caustic technique. In particular, the caustic mass profile 
assuming constant Tp somewhat overestimates the true 
mass profile within ~ 0.5/i _1 Mpc in simulated clusters 
( Serra et al.| 2011 1. We include these issues in our as- 
sessment o t the intrinsic unc ertainties and biases in the 
technique ( Serra et al.||2011 ). 

Some investigators have experimented with a modified 
caustic technique utilizing a radially-de pendent Tq (r) 
tailored to match simulated c lusters (e.g., Biviano & Gi- 
rardi|2003 Lemze et al.|2009 ) . Because one goal ot mea- 
surmg mass profiles with the caustics is to test the pre- 
dicted mass profiles from simulations, our approach is to 
assume a constant rather than impose a functional 
form for Tp(r). 

Note that the caustics extend to different radii for dif- 
ferent clusters. These differences result in part from the 
varying physical size subtended by the Hectospec field at 
different redshifts. 



D99 and Serra et al. (2011) show that the appearance 
of the caustics depends strongly on the line of sight; pro- 
jection effects can therefore account for most of the dif- 
ferences in profile shape in Figure [3] without invoking 
non-homology among clusters. 



We use the caustics to determine cluster membership. 
Here, the term "cluster member" refers to galaxies both 
in the virial region and in the infall region. Figure [2] 
shows that the caustics effectively separate cluster mem- 
bers from background and foreground galaxies. So me in- 



terlopers un avoidably lie within the caustics (e.g., [Serra 



et al. 2011). Serra et al. (2012, in preparation) applied 



the caustic technique to 3000 clusters extracted from N- 
body simulations. The technique identifies at least 95% 
of the true cluster members within 3r2oo- Only 2% of 
the galaxies inside the caustics and projected within r2oo 
are interlopers; the fraction of interlopers reaches 9% at 

3^200 ■ 

3.3. Comparison to Virial Mass Estimates 



Zwicky (1933 1937) first used the virial theorem to 
estimate the mass ot trie Coma cluster. With some mod- 
ificati ons, notably a corre ction term for the surface pres- 
sure (The & White 1986 ), the vi rial theorem remains in 
wide use (e.g., |Girardi et al. 1998 and references therein) . 
Jeans analysis incorporates the rad ial dependence of the 
projected velocity dispersion (e.g., Carlberg et al.| 1997 
|van der Marel et al.|[2000| |Biviano Jz Girardip003[ anc 
references therein) and obviates the need for a surface 
term. 

Virial mass estimators rely on the assumption that 
galaxy orbits are in equilibrium, an assumption that is 
certainly violated in the infall region and probably also 
in the inner regions of clusters with significant substruc- 
ture. Nevertheless, we apply the virial mass estimator to 
the HeCS clusters to check our caustic mass estimates. 
We must define a radius of virialization within which the 
galaxies are relaxed. We use r2oo as defined from the 
caustic mass profile (Table [3| and include only galaxies 
within the caustics. We thus assume that the caustics 
provide a good division between cluster galaxies and in- 
terlopers (see Figure J2j. 

We calculate the virial mass according to 



_ 3tt ajRpy 
2 G 



(3) 



where R PV = 2N(N - l)/Y,i,j>i R i^ is the projected 
virial radius and a 2 — J2i( v i ~ v) 2 / (N — !)■ If t ne system 
does not lie entirely within ^200, a surface pressure term 
3PV should be added to the usual virial theorem so that 
2T + U = 3PV. The virial mass is then an overestimate 
of the mass within r 2 oo by the fractional amount 



C = 47rr: 



P(^200) 



200 Pf-200 



f W0 Anripdr 



oy(r-2oo) 
o-(< r 20 o) 



(4) 
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024024024 

R p (h" 1 Mpc) 

Fig. 2. — Redshift versus projected radius for galaxies around HeCS clusters. The caustic pattern is evident as the trumpet-shaped 
regions with high density. The solid lines indicate our estimate of the location of the caustics in each clu ster. C lusters are ordered left-to- 
right and top-to-bottom by decreasing X-ray luminosity. This figure shows the first 12 clusters; Figures |16|19| show the remaining HcCS 
clusters. 



where oy(?*20o) is the radial velocity dispersion at r2oo 
and ct(< r 2 oo) is the enclosed total ve locity dispersion 
within r 200 ( e ;g-> ' 
ing cases of circu 



Girardi et aL| |1998| ). In the limit- 
ar, isotropic, and radial orbits, the 
maximum value of the term involving the velocity dis- 
persion is 0, 1/3, and 1 respectively. We estimate the 
uncertainties using the limiting fractional uncertainties 
7r _1 (21niV) 1 / 2 A r ~ 1 / 2 . These uncertainties do not include 
systematic uncertainties due to membership determina- 
tion. Table|3]lists the virial and caustic mass estimates at 
the radius r2oo determined from the caustic mass profile. 

Figure [4] compares the virial and caustic mass esti- 
mates at r2oo- The mean ratios of these estimates are 
M c /M v = 1.12 ± 0.04. The caustic mass estimates are 
slightly larger than virial mass estimates. Including a 
correction factor C « 0.1 — 0.2 M vir , consistent with the 
best-fit NFW profiles (see also jCarlberg et al ||1997[[GT 



rardi et al.||l99"8] |Koranyi fc Geller||2000| |Rines et a 



|2003| CIRS), would lead to a larger difference. For the 
G1RS clusters, this ratio is slightly below unity. When 
we combine the CIRS and HeCS cluster samples (130 
total), we obtain an average of M c /M v = l.Ollloioli- 
These results indicate that the caustic mass profile and 
the virial theorem yields similar masses on the scale of 
the virial radius (approximately r2oo). 

Figure [3] compares the mass profiles estimated from 
the caustics and the virial theorem. The virial mass pro- 
files are simply the virial mass estimator applied to all 
galaxies (inside the caustics) within that projected ra- 
dius. Because the virial theorem only applies in regions 
where the galaxies are in equilibrium, we only display the 
virial mass profiles in the range (0.75 — 1.3)r2oo (corre- 
sponding roughly to r 50 o to r 10 o)- The virial and caustic 
mass profiles generally agree. That is, near the virial 
radius, the caustic mass profiles do not appear to consis- 
tently overestimate or underestimate the mass relative to 
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Mpc) 



Fig. 3. — Caustic mass profiles for the HeCS clusters. The thick solid lines show the caustic mass profiles and the thin lines show the lcr 
uncertainties in the mass profiles. The inner vertical solid line in each panel shows the radius T200- The next vertical line shows the smaller 
°f r 5.6 (the limit of bound structure) and r max (the maximum radius where the caustics are detected). For clusters with r max < r'5.6, 
a dashed vertical line shows a lower limit on r§ .6 assuming no mass is present outside r max . Green shaded regions in dicate the virial 
mass profile in the range (0.75-f .3)r200 (approximately from rsoo to rioo). This figure shows the first 12 clusters; Figures |20|23| show the 
remaining HeCS clusters. 



the virial mass profiles. This result supports our use of 
caustic mass profiles as a tracer of the total cluster mass 
profile. 

3.4. Virial Masses and Ultimate Halo Masses 

The caustic mass profiles allow direct estimates of the 
virial and turnaround radius in each cluster. For the 
virial radius, we estimate r2oo- In our adopted cosmology, 
a system should be viri alized inside the slightly larger 
radius ~r wo k, 1.3r 2 oo (Eke et al. 19961. We use r2oo 



l-3r2oo ( Eke et al. 
because it is more commonly used in the literature and 
t hus allows easier compariso n of results. 
Nagamine & Loeb (|2003|) show that particles within 



a radius enclosing an overdensity of <5 C =17.6 (enclosed 

density 5.56p cr it or ^-p cr it) in the present epoch remain 
bound to the central halo in numerical simulations of the 



far-future evolution of large-scale structure in a ACDM 
universe. This criterion is furt her supported by simu- 

20051 



lations by other in vestigators dBusha et al. 2003 
Diinner et aL|2Tj06]) To be precise, IDiinner et al 



show that about 10% of the particles within this radius 
eventually become unbound, but that more distant par- 
ticles constituting 13% of the mass within <5 C =17.6 are 
ultimately accreted by the halo; thus, the mass within 
£ c =17.6 is 1.03 times smaller than the mass of the halo 
when the scale factor is a=10(r] If the w parameter in 
the equation of state of the dark energy (Pa = wpa) sat- 

6 In previous work, we used the slightly more generous defini- 
ti on of the turnaround ra dius rturn determined from equation (8) 
of |Regos fc Geller| fi989 1 assuming Qm = 0.3. For this value of 
!l m , the enclosed density is 3.5p c at the turnaround radius, versus 
5.58p c for <5 C =17.6. 
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10" 

IU. ft" 1 M 0> 

Fig. 4. — Caustic masses at r200 (determined from the caustic 
mass profile) compared to virial masses at the same radius. Solid 
squares show HeCS clusters and open squares show clusters from 
CIRS. Errorbars show 1<t uncertainties and the dashed line has 
slope unity. The solid line is the bisector of the ordinary least- 
squares fits. 



isfies w > — 1, the dar k energy has little effect on th e 
turnaround overdensity ( Gramann fc Suhhonenko|2002 1 . 
Varying fl m in the range (J. 1)2-1 only changes the in- 
ferred value of r turn by ±10%; the uncertainties in r turn 
from the uncert ainties in the mass p rofile are comparable 



or larger (D99; Rines et al 



quantify the ultimate mass o 



2002). Busha et al. 



(2005) 



dark matter haloes in their 



simulations by the ratio M5.6/M200, i.e., the ratio of the 
ultimate mass to the present value of M 2 qq. They find 
that the mass ratio follow a log-normal distribution with 
a peak at 2.2 and a dispersion of 0.38 (about 10% of halos 
occupy a high-end tail due to haloes merging with larger 
halos) . 

Table [3] lists r 2 oo, ^5.6, and the masses M200 and M5.6 
enclosedwithin these radii. For some clusters, the maxi- 
mum extent of the caustics r max is smaller than r 5 6 . For 
these clusters, M^.q and r$.Q are minimal values assuming 



that there is no additional mass outside 



The best 



estimate of the mass contained in infall regions clearly 
comes from those clusters for which r max > rs.6. The 
average mass within r^ ^ for these clusters is 1.99±0.11 
times the virial mass M 200 > in remarkable agree ment with 
the prediction of 2.2M 20 o ( |Busha et aT|l2005[ ). Further, 
the dispersion in In M max /M 2m is 0.30 (Figure 5|), similar 
to th e value of 0.38 for simulated clusters ( jBusha et al. 
2005 1 . The HeCS determination of M 5 . 6 /M 200 demon- 
strates that clusters are still forming in the present epoch 
(this ratio is larger than unity). Further, the measure- 
ment is consistent with the CIRS estimate of 2.19±0.18 
(the CIRS estimate refers to M3.5 rather than M 5 . 6 ; us- 
ing M5.6 for CIRS clusters would bring this value even 
closer to the HeCS estimate). 

The remarkable agreement of the HeCS estimate of 
a cluster's ultimate halo mass with the prediction from 
simulations is a new test of ACDM structure formation 
theory. 

Because estimating the ultimate halo mass of a clus- 
ter (assuming a ACDM model) requires determining the 




Fig. 5. — Ratio of mass M ma x within the maximum radius of 
the caustics (or r^ g) to the mass Af2oo within r2oo- Filled squares 
show clusters for which r max > r^ g and open squares show clusters 
with r max < ^5.6- The thick solid line shows the mean value of 
^5.6/^200 f° r clusters with r max > r5.fi. Thin solid lines show 
the uncertainty in this value. Triangles show CIRS clusters. The 
dashed line at Ms.6=2.2M200 is the ultimate mass of a halo in 
the far future (when the scale factor is q=10 0) compared to th e 
present-day mass M200 from the simulations of |Busha et al,| | |2005| ). 
The distribution of M(a = 100)/M2oo(c = 1) m their simulations 
is well-described by a log-normal distribution with a dispersion 
shown by the magenta hatching (sloping down to the right). The 
dispersion in ln(M5.g/M2oo) f° r the HeCS clusters is shown by the 
dense blue hatching (sloping up to the right). 



mass profile to a radius of r$^, the caustic technique is 
the only mass estimator that provides a direct probe of 
the ultimate halo mass. Weak lensing can detect shear 
signals at these radii, but the contribution of line-of-sight 
structure is difficult to separate from the lensing shear of 
the galaxies bound to the cluster. 

The agreement between the estimates of M5. q/M 2 qo 
and M3.5/M200 between the CIRS and HeCS samples 
suggests that the overall shapes of cluster mass profiles 
into their infall regions are not strongly dependent on 
cluster mass or redshift. 
In contrast, our mass profiles are per haps in tension 



with the analysis of |Tinker et al 
pares our estimates ot M max 



(2005) 



Figure Y6\ com- 
/M 2n n a s a fun ction of r max 



to the simulations of Tinker et al. (2005). Similar to 
CAIRNS and CIRS, the HeCS clusters tend to lie below 
the simulations, suggesting that either M200 is overes- 



ever, the ~50% offset suggested by |Tinker et al. (2005) 
is difficult to reconcile with the good agreement between 
observed and simulated values of M 5 , 6 /M 20 o discussed 
above (see further discussion in § 4). 

One striking result of this analysis is that the caustic 
pattern is often visible beyond the radius rs.6 that marks 
the maximum radius of galaxies that are gravitationally 
bound to the cluster. This result suggests that clusters 
may have strong dynamic effects on surrounding large- 
scale structure beyond the turnaround radius. For our 
assumed cosmology, most of this large-scale structure is 
not bound to the cluster. 
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TABLE 3 

HeCS Characteristic Radii and Masses 



Cluster 


»"500 


»"200 


7*5.6 


'f'm ax 




M 20 o 








M5.6 


M m ax /A/200 




Mpc 


Mpc 


Mpc 


Mpc 




10 14 M Q 




1O 14 M 




1O 14 M 








A267 


0.80 


1.19 


4.26 


4.85 


4 


95± 





31 


5 


20± 


0.44 


6 


32±0.44 


1 


28±0 


12 


Zwl478 


0.41 


0.64 


2.68 


3.64 





66± 





03 





74± 


0.10 


1 


38±0.08 


2 


09±0 


16 


A646 


0.67 


0.98 


4.24 


10.00 


2 


44± 





12 


1 


99± 


0.23 


5 


58±0.34 


2 


29±0 


21 


A655 


0.74 


1.08 


4.46 


8.89 


3 


33± 





16 


3 


47± 


0.32 


6 


50±0.42 


1 


95±0 


26 


A667 


0.73 


1.02 


4.09 


7.68 


2 


86± 





05 


1 


58± 


0.21 


5 


09±0.12 


1 


78±0 


06 


A689 


0.55 


0.80 


3.95 


5.66 


1 


55± 





05 





57± 


0.08 


5 


32±0.16 


3 


43±0 


21 


A697 


0.76 


1.13 


4.60 


3.74 


4 


42± 


2 


10 


3 


49± 
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3.5. Cluster Scaling Relations 

Scaling relations between simple cluster observables 
and masses provide insight into the nature of cluster as- 
sembly and the properties of various cluster components. 
Establishing these relations for local clusters is critical 
for future studies of clusters in the distant universe with 



the goal of constraini ng dark energy (Majumdar & Mohr 



2004 



Lin et al.||2004 |. 



We apply the prescription of Danese et al. ( 1980 ) to 



determine the mean redshift czq and projected velocity 
dispersion a p of each cluster from all galaxies within the 



caustics. We calculate a p using only the cluster members 
projected within r2oo estimated from the caustic mass 
profile. Note that our estimates of 7'2oo do not depend 
on o~p. 

Figure [7] shows the M 2 oo — <J P relation. The tight rela- 
tion indicates that the caustic masses are well correlated 
with velocity dispersion estimates. The good correla- 
tion is not surprising because both parameters depend 
on the galaxy velocity distribution. The best-fit slope is 
M200 oc (j2.90io.i5 with uncertainty estimated from 
jackknife resampling. The dashed line in Figure [7] shows 
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Fig. 6. — Ratio of mass M max within the maximum radius 
of the caustics (or rs.g) to M200 versus the maximum radius 
fmax of the caustics. Red hexagons are HeCS clusters with 
M200 > 3 X lO 14 /i _1 Af , red triangles are HeCS with M200 < 
3 X 10 14 /i — 1 Mq. Filled points have r max > rs 6- Black points are 
CIRS clusters. Filled squares show clusters for which Trnux > r 5.6 
and open squares show clusters with r max < rg.g. See Figurejfjjfor 
the typical uncertainties in Mmax IMm n. The lines show th e mass 
profiles predicted by the simulations of |Tinker et al.| ( |2005[ ) ■ 




1Q 13 JQ14 10 15 

M 200 ,,„, (h-' M Q ) 

Fig. 7. — Caustic masses at T200 compared to velocity disper- 
sions within T2oo- Squares and triangles show HeCS and CIRS 
clusters respectively. The solid line is the bisector of the ordinary 
least-squares fits. The dashed line is the a v — M2 00 relation of 
dark matter particles in cosmological simulations by |Evrard et al."| 
12008) ■ 



the M200 — ct p relation found by Evrard et al. 
dark matter particle s in simulated dark mat 
Evrard et alj (2008) find that this relation is 



([20081 for 
er haloes, 
insensi- 



tive to variations m cosmological parameters or numeri- 
cal resolution (above 10 3 tracer particles). The excellent 
agreement between the observed CIRS and HeCS clus- 
ters and the virial scaling relations from simulated dark 
matter haloes (slope 2.98±0.02) suggests that the caustic 
technique yields accurate mass estimates. 



Lx/E(z) (10" h- 2 erg S" 1 ) 

Fig. 8. — Rest-frame ROSAT X-ray luminosities compared to 
velocity dispersions within T200- Several outliers are labeled, and 
arrows indicate corrected luminosities for A689 and A1758. 

Figure [8] compares the HeCS velocity dispersions to 
rest-frame X-ray luminosities in the ROSAT band. There 
is a clear correlation, but there are also several outliers. 
Figure [S] labels several of these outliers. Two outliers, 
A689 and A1758, have problematic X-ray luminosities. 
Figure [8] shows that A689 lies close to the main locus 
of points when the X-ray luminosity is corrected for the 
central point source. We defer a full analysis of the Lx — 
a p relation to future work. 

3.6. Comments on Individual Clusters 

Clusters share many common features, but any large 
sample of clusters contains some complex systems. We 
comment on some of the most exceptional cases here. We 
present a comparison of caustic mass profiles and weak 
lensing mass profiles of several HeCS clusters in Gcllcr 
et al. (2012). 

• A667 The hierarchical center of A667 is located 
3Q7h- 1 kpc N of the BCG. The X-ray center is close 
to the position of the BCG. The spatial distribution 
of cluster members shows significant substructure 
to the North. This substructure accounts for the 
offset in the centers. 



• A689 This cluster was classified in BCS as a com- 
pact X-ray source, suggesting possible contamina - 
tion by a central point source ( |Ebelinget al.|1998[ ). 
A recent Chandra observation confirms that most 
of the X-ray luminosity in BCS is due to a central 
point source ident i fied a s a BL Lac (|Giles et al. 



2012[ ). |Giles et al. ( ]2012| ) estimate that the cluster 
luminosity in the BCS catalog is overestimated by 
about a factor of 10. With the revised luminosity, 
A689 lies below the flux limit of our flux-limited 
sample (Figure [T]). The Hectospec redshifts con- 
firm that the mass of A689 is significantly smaller 
than the masses of other clusters with similar (un- 
corrected) Lx- Figure [8] shows that A689 is not an 
outlier in the Lx — o~ p diagram when using the cor- 
rected Lx- This cluster highlights the importance 
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Rines et al. 




Fig. 9. — Top panels: Redshift versus projected radius for 
MS0906 (left) and A750 (right). The caustic pattern of A750 is 
much clearer when the plot is made with A750 at the center (cf. 
Figured. Bottom panels: Mass profiles of MS0906 (left) and A750 
(right) LJ 

of high-resolution X-ray observations in determin- 
ing accurate X-ray luminosities. 

• MS0906/A 750 This pair of clusters is a curious sys- 
tem. MS0906+H is an X-ray cluster at z=0.1767 
detected in the Einstein Me dium Sensitivity Slew 

A750 is a cluster 



1992) 



Survey (Henry et al 
at z=0.1640 located only - !)' (0.63/i _1 Mpc) away 
fr om the X-ray center of MS09 06 (see Figu r e 3.39 
of iMaughan et al.|]2008| ). |Carlberg etat] (| 1996 1 
noted that MS0906 "appears to be an indistinct 
binary in redshift space". With the denser sam- 
pling of HeCS, the infall patterns of the two clus- 
ters are separable. Figure [9] shows the caustics 
and mass profiles of these two clusters. The two 
clusters are separated by 3250km s _1 (rest-frame), 
suggesting that they are not gravitati onally bound. 
Note th at the weak lensing map of |Okabe et al.| 
(2010a) shows two distinct mass components cen- 
tered approximately on MS0906 (component A750- 
C in their Figure 32) and A750 (component A750- 
NW1 in their Figure 32); MS0906 has larger sur- 
face mass density and A750 has a larger luminosity 
density of red-sequence galaxies. The X-ray lumi- 
nosity of MS0906 is mu ch larger than that of A750 
(Maughan et al. 2008). Our caustic mass profiles 
indicate that the two clusters have roughly equal 
mass. This complex system indicates that the re- 
lation between cluster mass and different observ- 
ables is complicated. A more detailed analysis of 
this system is presented in Geller et al. (2012). 



A 773 A773 contains a radio halo (jGiovannini et al." 



1999), a f eature commonl y assoc iated with major 



mergers. [Barrena et al. (|2007|) studied the dy- 



namics of 100 cluster members and found evidence 
of complicated dynamics, including two velocity 
peaks separated by « 2500km s~ . Each velocity 
peak contains one of the two bright galaxies ly- 



ing in the cluster center. A third galaxy (SDSS 
J091758.60+515104.6) is intermediate in bright- 
ness to the two bright galaxies in the core; this 
third galaxy is in the lower-velocity peak and is 
relatively isolated. The caustics contain both sys- 
tems, although some galaxies in the high-velocity 
peak lie outside the caustics. Because the caustic 
method does not require equilibrium, the caustic 
method may be a more robust mass estimator for 
A773 than virial analysis. The hierarchical cen- 
ter has a redshift similar to the low-velocity peak, 
and the caustics at large radii are centered at the 
same redshift. Note that our estimate of the veloc- 
ity dispersion (lllO^Qkm s _1 ) is smaller than the 



-70 r 

global velocity dispersion 1394lggkm s~ 



of 



Bar- 



rena et al. ( 2007 ) , although they find much smaller 
velocity dispersions of individual subclusters. 



A1437 This cluster was studied by |Pimbblet et al.| 
(2006). They measured a velocity dispersion of 
1152l^km s _1 , smaller than but consistent with 
our estimate (a p = 1233igi 2 km s _1 ). 

A 1689 T his cluster is a famous stro ng -lensing clus 



ter (e .g., Broadhurst et al. 2005). 



Lokas et al. 



(2006) studied the dynamics of Al689 from liter 



ature data and find significant line-of-sight sub- 
structu re. The dynamics h ave also been stud- 
ied by |Lemze et al. (20091 using more exten- 



sive V L T/ V1MUS spectros copy obtained by Czoske 
( 2004 ) . Lemze et al. ( 2009 ) show that the mass pro- 
file determined with the caustic technique is con- 
sistent with mass profiles from gravi tational tens- 
i ng, X -ray data, and Jeans' analysis. |Haines et al.| 
( |2010[ ) use 1009 Hectospec redshifts (from an in- 
dependent investigation) to study the properties of 
star- forming galaxies d etected by Herschel. Note 
that |Lemze et al. ( 2009 ) claim to detect the "edge" 
of A 1689 by noting a sharp decrease in galaxy den- 
sity at R p 2.1/i -1 Mpc. Figure ^1 shows sev- 
eral spectroscopically confirmed members beyond 
this radius, and we detect caustics extending to 
R p = SJli^Mpc. 

• A 1758 This system is a merger of multiple X- 
ray clusters at z=0.28. Several authors se parate 
the cluster into A1758N and A1758S (e.g ., lEbeb] 



ing et al. 1998 David & Kempner| |2004|) . The 



combined flux of the two clusters (jBohringer et al. 



2000b ) exceeds our flux limit; the individual fluxes 
lie below the flux limit. David & Kempner (20041 



showed that A1758N is itself a merger of two 7 keV 
clusters and that A1758S shows evidence of a re- 



cent merger. Okabe k, Umetsu| ( |2008[ ) used Subaru 
to analyze the lensing properties of A1758N and 
A1758S and found confirming evidence that both 
components are undergoing mergers. Ra gozzme 



et al. (20121) use a higher- resolution lensing anal- 



ysis of A1758N and conclude that A1758N con- 
sists of two separate clusters A1758N:NW and 
A1758N:SE. Because A1758 is one of the highest- 
redshift clusters in HeCS, it is not sampled very 
deeply. A detailed dynamical analysis of this com- 
plex system would require additional data. The 
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complex dynamics of the mergers in A1758 may 
contribute to its unu s ually high star formati on rate 



(Haines et al. 2009). 



(20121 study 



Boschin et al. 
the dynamics of A1"75BN:NW and Al758rT:SE us- 
ing TNG spectroscopy. They find a much larger 
velocity dispersion (a p ~ 1300km s _1 ) than we do, 
but they note that many high- velocity-offset galax- 
ies might be members of subclusters. 

A 1902 The hierarchical center is coincident with a 
close grouping of galaxies 531/i _1 kpc E of the BCG. 
The close grouping is probably a merging event, 
leading to a very high binding energy that causes 
the hierarchical analysis to choose this grouping 
as the cluster center. The X-ray center is located 
close to the BCG, and the spatial distribution of 
cluster members shows many more members East 
of the BCG than West of the BCG. The BCG has 
some close companion galaxies that may have been 
missed by either SDSS photometry or by fiber colli- 
sions. In cither case, the BCG companion galaxies 
would be omitted from the redshift catalog used for 
the hierarchical analysis; including some of these 
companions could conceivably relocate the hierar- 
chical center close to the BCG. 

A 1914 Although the X-ray contours of A1914 are 
fairly regular, it contains a large radio halo and 
significant X-ray temperature anisotr opy, two fea- 
tures consistent with a major merger ( Govoni et al. 
2004 ) . Similarly , a weak lensing map by |(Jkabe 
&: Umetsu ( 2008[ ) shows significant substructure m 
both the lensing mass distribution and the luminos- 
ity distribution (which includes two bright galaxies 
close to but not coincident with the X-ray center). 
The HeCS data (Figure [2]) show a clear caustic pat- 
tern but with an unusual shape: the median red- 
shift of cluster members increases from the center 
to ~ 0.5/i _1 Mpc, then decreases from this radius 
to ~ 2.0/i~ 1 Mpc. Inspection of the spatial distri- 
bution of the cluster members reveals that a large 
group is located ps 1.5/i -1 Mpc NW of the X-ray 
center. The group is offset by s» —1000km s _1 
from the redshift center of A1914 and has a ve- 
locity dispersion of a p ~ 350km s _1 . This group 
largely explains the unusual appearance of the in- 
fall pattern in Figure [2j We defer a full dynamical 
analysis of this complex system to future work. 

A 1930 The hierarchical center is located 
385/i _1 kpc from the X-ray center. The hier- 
archical center is 3f 5 SSE of the BCG and V from 
the second-ranked galaxy. The X-ray center is 2'5 
WSW of the BCG. The BCG is located about 1.'5 
NNE of the midpoint of the hierarchical center 
and the X-ray center. 

RXJ1504-03 This cluster is near the edge of the 
DR6 photometric footprint. Our target photom- 
etry is based on DR5, which did not contain the 
entire cluster. Part of the Hectospec field was con- 
tained in DR5 imaging that did not yield photom- 
etry. We used the SDSS Navigate tool to identify 
and include red and blue galaxies in our target list 
by inspecting the unprocessed SDSS imaging. A 



small portion of the Hectospec field was covered by 
neither SDSS photometry nor Atlas imaging. We 
added targets from POSSII plate scans in this re- 
gion. 

RXJ1504 is well know n as a cooling-core cluster 
(Bohringer et al. 2005), and it is the most X-ray 
luminous cluster in HeCS. Vikhlinin et al (2009) 
estimate a large mass for RXJ1504 based on Chan- 
dra observations and the Yx estimator. However, 
Figures [2] and [8] show that its velocity dispersion 
is smaller than many clusters of comparable X-ray 
luminosity. This result suggests that X-ray mass 
estimates of RXJ1504 (including Yx) may be bi- 
ased due to the cooling core. Recently, [Zhang 



et al. (2012) studied the dynamics of galaxies in 
RXJi50¥"using VLT/VIMOS spectroscopy. They 
identify 53 cluster members and compute a ve- 
locity dispersion of <7 p =(1132±95)km s -1 , a value 
significantly larger than ours (<7 p =779^75 5 km s _1 , 
Table [I]). Comparing their Figure 7 to our Fig- 
ure [2] indicates that their membership selection is 
more generous than our caustic selection; HeCS in- 
cludes more redshifts (120 members versus 53) and 
covers a wider field-of-view t han the spect ra con 
tained in|Zhang et al. (2012). In addition, Zhang 



et al. (2013?) identity several blue galaxies as clus 
ter members: if these galaxies have a larger ve- 
locity dispersion than the red galaxies we target, 
this difference could partly explain the difference 
in measured velocity dispersions. However, anal- 
ysis of blue galaxies in other HeCS clusters ({ 5.3 1 
shows no significant color dependence of the ve- 
locity distribution of members identified with large 
redshift samples and the caustic technique. 



• A2055 This cluster was studied by |Pimbblet et al.| 
(2006). They measured a velocity dispersion of 
1046lggkm s _1 , significantly larger than our es- 
timate (<7 p = 676lg4). Inspection of the infall 
pattern shown in their Figure 5 reveals that they 
classify three galaxies with large velocity offsets 
as members; similar galaxies are classified as non- 
members by the caustic technique (Figure [2]) . 

• A2219 The kinematics of this well-know n lensing 
cluster ( Smail et al.|1995 1 were studied by Boschin 
et al I pU4| " using 132 redshifts within 5' of the 
BUG. They find significant evidence of substruc- 
ture and possible merging activity. We include 
their redshifts in our caustic analysis and we re- 
moved these galaxies from our Hectospec target 
list. A 2219 contains a radio halo (Giovannini et al. 
|1999[ ). Figure [2] shows that there are a few galax- 
les projected in front of A2 219; while we classify 



Boschin et aL| ( |2004 1 
'his membership dir 



these galaxies as foreground 
classify them as members 
ference probably accounts for most of the differ- 
ence between their estimated velocity dispersion 
(cr p = 1438lgg 9 km s _1 ) and our estimate (Table 
[T|). These foreground galaxies could also enhance 
the probability of observing strongly lensed back- 
ground galaxies. 
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• Zw8197 The hierarchical center is located 8' (667 
kpc) W of the X-ray center. The X-ray center lies 
close to the BCG, but there are few/no galaxies 
around the BCG. The spatial distribution of clus- 
ter members is remarkably flat; the BCG lies in a 
sparse region of cluster members. The large magni- 
tude ga p between the BCG a nd other cluster mem- 
bers led Santos et al. (2007) to identify Zw8197 as 
a candidate fossil group. 

• A2261 This system is near the edge of the DR6 
photometric footprint. Although imaging is avail- 
able through the SDSS Image List tool, photome- 
try was not available in DR6 for much of the clus- 
ter (this situation can arise when a field is taken 
in poor seeing). We used the Guide Star Catalog 
to identify likely galaxies in the cluster region and 
the Image List tool to visually inspect these candi- 
dates. We experimented with using colors from the 
Guide Star Catalog photometry, but these identi- 
fications are less reliable than visual classification 
of the SDSS thumbnails. DR8 contains photome- 
try for A2261. We used this photometry to select 
targets for an additional Hectospec pointing. We 
increased the sampling density and included blue 
galaxies to examine detailed issues in the strong 
and weak lensi ng determination of the cluster mass 
( |Coe et all[2012| . 

The D99 binary tree analysis we use for all other 
clusters locates the center of A2261 on a struc- 
ture 613 (« l/i^Mpc) and 400km s" 1 away from 
the BCG. This center lies atop a tight grouping 
(30") of four bright red cluster members (two with 
L ~ 10-L*). However, the slightly different algo- 
rithm for cu tting the binary tree described in |Serra| 
et al. (2011) yields the cluster center on the BCG. 



These two different results are a consequence of the 
complex dynamics of A2261. We adopt the center 
closest to the BCG. The caustic mass profile is in- 
sensitive to the adopted center. 

Two teams have investigated the mass profile of 
A2261 using gravitational lensing. Okabe et al. 
p0T0a| find that A2261 is "over-concentrated" rel- 



ative to the expec tations from num erical simula- 
tions. By contrast, Coe et al. (20121 use HST lens- 
ing data to measure a concentration cl oser to the 
expect ed value for a massive cluster. |Coe et aL] 
(20121 show that lensing estimates of M200 can 
differ by ~25% depending on the assumed clus- 
ter geometry (spherical versus triaxial). This flex- 
ibility allows the lensing mass to agree with vari- 
ous X-ray hydrostatic mass estimates that differ by 
35% (and could be affected by non-thermal pres- 
sure support). However, both the X-ray and lens- 
ing mass estimates are larger than our estimate of 
M 200 (see Figure 12 of |Coe et al."||2012| ) . 



4. ENSEMBLE CLUSTERS: THE CLUSTER MASS 
PROFILE 

The largest uncertainty in determining dynamical 
masses of clusters is the influence of projection effects. 
Stacking clusters together to create an ens emble clus- 
ter can significantly reduce this uncertainty ( |Serra et al. 



2011 1. We create two sets of ensemble clusters using 
two methods of grouping the clusters. For all ensemble 
clusters, we eliminate A1758 and MS0906/A750 because 
they are double clusters (see ^3.6). 

First, we divide the clusters into quartiles of M 2 oo as 
determined from the caustic mass profile. We then stack 
the clusters in physical units (positions in kpc, velocities 
in kms^ 1 ). This method has the disadvantage that it 
relies on the parameters from the caustic mass profiles 
to assign the clusters into quartiles, so the resulting en- 
semble properties are not completely independent of the 
caustic mass profiles of the individual clusters. 

Second, we divide the clusters into quartiles of X-ray 
luminosity from the original ROSAT catalogs. We then 
stack the clusters in physical units (positions in kpc, ve- 
locities in fcms -1 ). This approach avoids any use of the 
individual cluster caustic mass profiles in the stacking 
parameters. 

Figure[l0]shows the two sets of ensemble clusters. Both 
sets show a clear distinction between cluster members 
and field gal axies in all quartiles. 

Figure [IT] shows the density profiles of the ensemble 
clusters. The caustics trace the mass profiles across 
nearly four orders of magnitude in density. The blue 
dashed line shows an NFW profile with a concentration 
of c = 3.6. The low-concentration NFW profile is clearly 
a much better fit to the inner region of the ensemble clus- 
ter. Concentrations of c «3-5 are predicted by num erical 
simulation s of cluster-sized dark matter haloe s (e.g.JBul- 



lock et al.||2001[ ). Further, [Serra et al.| ( f2011[ ) show tnat 



cluster mass profiles in cosmological simulations follow 
the extrapolation of an NFW fit (performed on the in- 
ner l/i _1 Mpc using the caustic technique) out to 3-4r 2 oo- 
That is, the simulated clusters follow NFW profile s well 
into their infall regions (see also Tavio et al. 112008 1. 



Some investigators have suggested that the high con- 
centration parameters found in cluster lensin g profiles in- 
dicate probl ems with ACDM cosmology (e . g . , | B roadhur st I 
To ameliorate projection effects, |Umetsu| 



et al. 



et al 



2008). 
WT) 



constructed a stacked cluster and deter 
mined the mass profile from both strong and weak grav- 
itational lensing. The magenta dashed line in Figure 11 
shows an NFW profile with concentration c = r2oo l r s = 
7.68; this pro file is an excellent fit to the stacked lens- 
ing cluster of Umetsu et al. (2011), one of the highest 



precision lensing profile estimates. Figure [TP] also shows 
the best-fit N FW profile (c=3.58) f rom the stacked lens- 
ing cluster of Okabe et al. (2010b) created from X-ray- 
selected clusters! The density profiles of the HeCS en- 
semble clusters generally lie between the lensing profile 
from the X-ray selected clusters and the profile from 
strong-lensing selected clusters. 

Our results for the ensemble HeCS clusters indicate 
little tension between the simulations and mass profiles 
determined from the caustic technique. One possible 
concern in measuring c for an ensemble cluster is that 
mis-centering could artificially smooth out the central 
density peak and thus artificially lower the measured con- 
centration. However, the offset between the hierarchical 
centers we use and the X-ray ce nter s is almost always 
less than 300/i _1 kpc, and Figure 11 shows a significant 
difference even at 500/i _1 kpc. Counteracting this effect, 
our assumption of constant overestimates the central 
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masses of simulated clusters ( Serra et al.|2011 1. We will 
investigate possible systematic uncertainties m measur- 
ing concentrations with caustics in future work (Diaferio 
et al. in preparation). 

The lower panels of Figure [TT] show the cumulative 
density profiles of the ensembleclusters. These profiles 
determine the values of the radii r$ , where the enclosed 
density within rg is equal to Sp c . For all quartiles, the 
caustic mass profiles extend well beyond r2oo, but not 
quite to r^.e (the maximum radius of particles bound to 
the cluster in the far future). This limitation is primar- 
ily due to the limited radial extent of the HeCS data: 
a radius of 30' corresponds to 4.2/i _1 Mpc at z = 0.2; 
the minimum value of r^^ for the highest-Lx ensemble 
cluster is 5.8/i _1 Mpc. 

Figure [12] shows the residuals from NFW fits to the in- 
ner l/i _1 Mpc of the ensem ble clusters ( fits pe rformed on 



Figure [14] shows the absolute magnitudes of cluster 
members versus their redshifts. The solid line show 
M r = —20.60 + 51ogh, approximately the characteristic 
magnitude M* of the luminosity fu nction of field galax- 



the caustic mass profile). Serra et al. (2011) show that 



this procedure yields accurate predictions of their simu- 
lated cluster mass profiles. The density profiles are noisy 
but seem to indicate that observed densities at large radii 
are slightly smaller than the densities of the extrapolated 
NFW fits. The lower panels of Figure [12] show that the 
cumulative density profiles agree within ~20% of the ex- 
trapolated NFW profiles far beyond T2oo- The smallest- 
Lx quartile and the two smallest-M2oo quartiles have 
smaller cumulative densities than the extrapolated NFW 
profiles. We speculate that this deficit is a combination of 
the low-c values of these profile fits (c ~2-3) and a dearth 
of observed galaxies at large radii (Figure 10 ) resulting 
from the lower average redshifts of the constituent clus- 
ters and the finite field-of-view of Hectospec. Overall, 
the agreement between caustic mass profiles and NFW 
profiles extrapolated to large radius is excellent. This 
dynamical agreement is independent support of similar 
results deriv ed from weak lensing profiles extending to 
large radius ( Okabe et al.||2010"a Umetsu et al.||2011 ). 



5. PROPERTIES OF CLUSTER GALAXIES 

One goal of HeCS is to study the spectroscopic proper- 
ties of galaxies in clusters and their infall regions. We de- 
fer a complete analysis to future work. Here, we investi- 
gate the possible impact of our target selection algorithm 
(which favors red-sequence galaxies) on the estimates of 
cluster mass profiles. 

5.1. Red Sequence Galaxies 

Because we targeted galaxies on the red sequence (Ap- 
pendix), many o f th e cluster members lie on the red se- 
quence. Figure 13 shows color-magnitude diagrams of 
cluster members. We computed K-corrections us i ng th e 
purely empirical K-corrections of Westra et al. (2010). 
One striking feature of Figure [13] Is the small scatter in 
the red sequence. The outer lines are offset from the red 
sequence by ±0.3 mag, approximately the color range 
where we assigned highest targeting priority. The inner 
lines are offset by ±0.1 mag, demonstrating that most 
cluster members lie within this much narrower range of 
color. This result suggests that the galaxy properties are 
very similar and that the SDSS photometric uncertain- 
ties are minimal. Lines with a slope of -0.04 (observed 
colors) or -0.025 (rest-frame colors) provide a remarkably 
accurate description of the red sequences. 



ics in SDSS (Blanton et al.||2003b|). The dashed line is 
rde 



one magnitude fainter [M* ± 1), a lumin osity limit of- 
ten a dopted for describing galaxies (e.g., Tinker et al. 
2005 1 . CIRS shows that this limit is the minimum depth 
for a cluster sample that yields enough members to iden- 
tify caustics. Figure [Ti] shows that all HeCS clusters at 
z < 0.25 meet this condition. The remaining clusters are 
still well-sampled because these clusters have larger X- 
ray luminosities (and hence more bright member galax- 
ies) than the CIRS clusters (Figure [IF. 

5.2. Extremely Red Cluster Galaxies 

Figure [15] shows that very few cluster members have 
g — r colors significantly redder than the red sequence. 
In fact, all of the six cluster members in RXJ2129 and 
A2261 with colors redder than the red-sequence cut have 
incorrect colors due to their proximity to nearby bright 
stars or galaxies. Using fiber magnitudes, all six have 
g — r colors below the red-sequence cutoff. 

Because the red sequence represents some of the old- 
est known stellar populations, colors redder than the red 
sequence would most likely arise from extreme dust red- 
dening. Thus, the absence of cluster members with colors 
more than 0.3 mag redder than the red sequence indicates 
that very few cluster members are extremely reddened. 

5.3. Blue Cluster Galaxies: Are They Common? Do 
They Affect Dynamical Mass Estimates? 

To test the sensitivity of our dynamical estimates on 
the red-sequence target selection, we observed several 
Hectospec pointings in the clusters A267, A2261, and 
RXJ2129 (all at z«0.23) to sa mple blue galaxies The 
Hectospec survey of RXJ1720 (Owers et al. 2011) also 
samples both red-sequence and blue galaxies. Figure 15 
shows that very few blue galaxies ar e cluster members , 
consistent with previous studies (e.g., |Rines et al.f 2005). 
The solid lines show the approximate limits ot the HeCS 
priority target selection used for the other HeCS clusters. 
The dashed line slightly bluer indicates our limit for low- 
priority targets. Very few galaxies below the solid lines 
(filled blue squares in the left panels) are cluster mem- 
bers. Combining the four clusters, 136 of 1050 cluster 
members (13.0%) lie blueward of our red-sequence cut. 
A two-sample K-S test indicates that the velocity dis- 
tributions of red and blue galaxies are consistent with 
being drawn from the same parent population. Further, 
the velocity dispersion of the ensemble cluster (all galax- 
ies) is only 0.7% larger than the velocity dispersion of 
the red-sequence galaxies. Limiting the sample to galax- 
ies inside r 2 oo, 37 of 487 members (7.6%) lie blueward of 
the red-sequence cut; including non-red-sequence galax- 
ies increases the velocity dispersion by 0.3%. We thus 
conclude that targeting red-sequence galaxies produces 
no significant bias in our estimates of dynamical masses 
or velocity dispersions. Previous claims of velocity segre- 
gation often relied on much smaller samples where mem- 
b ership classification m ay have been less robust 



Mahajan et al. (2011) recently studied the dependence 
of cluster galaxy properties on their projected velocity 
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Fig. 10. — Redshift versus projected radius in ensemble clusters. (Left) Quartiles of Lx (Right) Quartilcs of Af200- Quartiles are shown 
from top to bottom by decreasing mass/luminosity. Solid lines indicate the positions of the caustics for the ensemble clusters. 



relative to the cluster center. They bin galaxies both 
by projected radius and by velocity offset (e.g., galaxies 
with |Au|=(0-l)(Tp, l-2a p , 2-3<r p ). At fixed projected ra- 
dius, their Figure 3 shows little dependence of the color 
distribution on velocity offset, suggesting that the veloc- 
ity distribution does not depend dramatically on galaxy 
color, consistent with our results. 

Our Hectospec (and SDSS) redshifts show that some 
cluster members are "blue cloud" galaxies with colors in- 
dicating recent star formation (Figure 15). The fraction 
of "blue cloud" galaxies seems to increase with incr easing 
absolute magnitu de, similar to field galaxies (e.g. 



ton et al. 2003a). This conclusion could be tested with 



Blan- 



more extensive spectroscopic sampling of faint galaxies 
that lie blueward of the HeCS red sequence cuts. 

6. CONCLUSIONS 

We present 21,314 redshifts from the Hectospec Clus- 
ter Survey (HeCS). HeCS is a MMT/Hectospec spec- 



troscopic survey of X-ray-selected clusters contained in 
the imaging footprint of SDSS DR6. Our redshifts con- 
firm that infall patterns known as "caustics" are clearly 
present in X-ray clusters at moderate redshift. Combined 
with SDSS data, we define a sample of 10,275 cluster 
members. 

We use the infall patterns to compute mass profiles for 
the clusters extending in many cases to the turnaround 
radius of the cluster. In numerical simulations of a 
ACDM universe, the mass within the radius r^ ^ (the 
average density inside r^^ is 5.6p c (z)) is approximately 
equal to the ultimate mass of the cluster at late times. 
These simulations predict that this ultimate halo mass is 
w2.2M 2 oo- The HeCS mass profiles provide an observa- 
tional estimate of M5.6 = (1.99 ± 0.11)A/2oo, in excellent 
agreement with the predictions. 

The caustic technique enables a unique measure of the 
large-scale behavior of cluster mass profiles. The striking 
agreement with the theoretical predictions ( Nagamine & 
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Fig. 11. — Density profiles of ensemble clusters. Black, red, green, and blue solid lines show decreasing quartiles of Lx (left panels) 
or M'20t) (right panels). Top panels The magenta dashed line shows the JNFW prohlcs that best ht the stacked lensmg cluster ol Umetsu 
ct al. (2U11 strong-lcnsing selected). The blue dashed line shows an 1N1VV prohlc with concentration e=3.6, the best-ht prohle from the 
stacked lensing cluster ot Ukabc ct al. (2Ul(Jb X-ray selected). The lensing prohlc Irom the X-ray selected clusters better matches the 
caustic density prohlcs than the prohlc Irom strong- lensmg selected clusters, consistent with strong lensmg bias (Broadhurst et al. 2UU8). 
Horizontal lines indicate the critical density p c and p = U.'ip c - The density prohlcs detect signal out to relatively small densities. Bottom 
panels Cumulative density prohlcs of the ensemble clusters. Horizontal lines indicate enclosed densities of (500,200,5.6),Oc, where p c is 



evaluated at 2=U.lt>, the median rcdshitt ot the HcCJS sample. 



Loeb||2003| |Busha et aL"1|2005| [Dunner et al.||200 6) is an 

interesting new test ol models tor the growth ol structure 
in the ACDM cosmology. 

The ensemble clusters we construct here average over 
projection effects. The density profiles of these ensem- 
ble clusters closely resemble NFW profiles out to radii of 
at least 2r2oo ■ The NFW concentrations of our ensemble 
clusters are close to the predictions of simulations and do 
not appear to confirm the "over -concentration" problem 
of strong-lensing clusters (e.g., |Broadhurst fc Barkana| 
2008 1 . Future work will compare these caustic mass pro- 
files to those predicted by numerical simulations of dark 
matter haloes (Diaferio et al. in prep). 

The CAIRNS and CIRS projects demonstrated that 



caustic patterns are present in nearly all rich, X-ray lu- 
minous galaxy clusters. Here we show that these patterns 
are also prominent in clusters at z=0. 1-0.3. The ensem- 
ble clusters we construct enhance the visibility of these 
patterns. 

The HeCS clusters show a relation between Lx and a p 
(we will investigate the specific relation in future work); 
we show that outliers in this relation can be attributed to 
contamination from X-ray point sources (A689) or strong 
cooling cores (RXJ1504). The scaling relation M200 — o v 
between virial masses and line-of-sight velocity disper- 
sions is in excellent agreement with the scali ng relation of 
dark matter particles in simulated clusters (Evrard et al. 
2008 1. Other projects planned with HeCS are investi- 
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Fig. 12. — Residuals of NFW fits to the density profiles of ensemble clusters from Figure |11| Black, red, green, and blue solid lines 
show decreasing quartiles of Lx (left panels) or M200 (right panels). Top panels Residuals of trie ensemble density profiles from NFW fits. 
Note that the NFW fits are performed with the mass profiles within 1 Mpc and not to the density profiles. Bottom panels Residuals 
of cumulative density profiles of the ensemble clusters. NFW fits are performed using the caustic mass profile within lh~ x Mpc and 
extrapolated to large radii (see |Serra et al,|2011| ). 



gations of cluster scaling relations between galaxy dy- 
na mics and other ma ss probes (richness, Yx, Ysz,Miens, 
Rines et al.|2010 1 , a determination of the mass func 



see 



tion of clusters (eTg^ |Rines et al.||2008 



2009||Mantz et al.|2008||Henry et al.|20U9[ ) and a detailed 



Vikhlinin et al. 



study ot the photometric and spectroscopic properties of 
cluster galaxies. A companion paper (Geller et al. 2012) 
compares caustic mass profiles to those determined from 
weak gravitational lensing, the only other mass estimator 
that applies to the non-virialized infall regions surround- 
ing clusters. 
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Fig. 13. — Color-magnitude diagram of member galaxies of the 
HeCS clusters including K-corrections. Long-dashed lines show 0.3 
mag away from the red sequence, approximately the limits of our 
priority target selection. Most of the HeCS members are within 0.1 
mag of the red sequence (short-dashed lines), indicating that the 
color selection includes the vast majority of red-sequence cluster 
galaxies. The dotted line shows the color range 0.2 mag blueward of 
the red-sequence cut used to select additional targets to fill unused 
fibers. 



Fig. 14. — Absolute magnitude versus clustroccntric radius for 
cluster members in the HeCS clusters. Solid and dashed lines in- 
dicate absolute magnitudes of M* and M* + 1 respectively. 
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APPENDIX 
SPECTROSCOPIC TARGET SELECTION 

We selected targets for spectroscopy based on photometry from SDSS DR6. First, we extracted galaxy catalogs 
within 32' from the SDSS SQL serveiFl For each cluster, we identify the red sequence from a color-magnitude diagram 
using g — r colors and r-band apparent magnitudes. Following guidance from the SDSS web page^_J we adopt the 
composite model magnitudes for our apparent magnitudes and colors. Composite model magnitudes are a linear 
combination of the best-fit exponential and deVaucouleurs model profiles. Composite model magnitudes should have 
higher signal-to-noise than Petrosian magnitudes, especially for relatively faint galaxies such as those studied here. 

To identify the red sequence, we used a fixed slope of -0.04 magmag -1 for all clusters and chose the intercept 
based on visual inspection of the color-magnitude diagrarrj^] To eliminate some stars and artifacts (e.g., portions of 
diffraction spikes from bright stars), we restrict the sample to galaxies with r-band surface brightness brighter than 
22.9 magarcsec -2 (we later inspect objects with fainter surface brightnesses to include some galaxies eliminated by 
this cut). 

The Hectospec fiber assignment software xfitfib^j allows the user to assign rankings to targets. If we ignored the 
spatial positions of our targets, fiber collisions would prevent many objects in the central portion of the field (usually 
the center of the cluster) from being observed. We therefore assigned highest priority to galaxies brighter than a cluster- 
dependent limiting magnitude within ±0.3 magnitudes of the red sequence and within 2.4' of the X-ray center. We give 
second priority to red galaxies near the center and up to one magnitude fainter. Priorities are then assigned in annuli 
(outer radii 7'.5 and 15') according to apparent magnitude. In general, our target catalogs contained approximately 
twice as many targets as the number of fibers available. Most bright targets in the inner 15' are assigned fibers, as 
are many fainter targets in this region. Unassigned targets are typically in the outer parts of the Hectospec field, with 
a modest bias against targets in local overdensities due to fiber collisions. Later analysis shows that the ±0.3 mag 



7 http:/ /cas. sdss.org/astrodr6/en/tools/search/sql. asp 

8 http:/ /www. sdss.org/dr6/algorithms/photometry.html 

9 For 8 clusters observed in 2007 (A2111, A2187, A2219, A2259, 
A2396, A2631, A2645, and RXJ2129), we used a simple g-r color 



cut rather than a red-sequence cut. This small sampling difference 
has no significant impact on our results. 

10 https : //www . cf a. harvard . edu/- john/xf itf ibs/ 
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band around the red-sequence is significantly larger than the actual thickness of the red-sequence (see ( 5.1 1; thus, the 
precise choice of red-sequence intercept has essentially no impact on our results. Because cluster members are centrally 
concentrated, we found that prioritizing central objects is critical to obtain reasonably uniform sampling as a function 
of radius. 

Because fiber placement constraints would sometimes prevent fibers from being deployed, we add galaxies up to 0.2 
magnitude bluer than our red-sequence cutoff to the target list. We use the SDSS Image List tooP_Hto visually inspect 
all targets prior to observation. This visual inspection identifies artifacts and bright stars. We added a small number 
of targets that are not identified as separate objects by the DR6 photometric pipeline (these targets are usually merged 
with a nearby star or galaxy). 

We observed two Hectospec pointings per cluster for clusters at z>0.15 and one pointing per cluster for clusters at 
z=0.10 — 0.15. We adopted this strategy because SDSS redshifts of Main Sample galaxies at r<17.77 provide many 
members for the lower-redshift clusters (although not enough to enable a full caustic analysis). Also, the lower-redshift 
clusters have generally smaller X-ray luminosities (because HeCS is flux-limited) and thus have smaller expected 
masses. Combined with the larger angular size of these clusters, one Hectospec pointing combined with SDSS redshifts 
is generally sufficient to avoid sampling bias. For clusters at z>0.15, we chose apparent magnitude limits for each 
cluster to select 900-1300 candidate targets for two Hectospec pointings (potentially up to 600 targets). The bright 
magnitude limit for these clusters is 18. 5< r <20.0, corresponding to M* + 1 or fainter for most clusters (Figure 



14). For clusters at z =0.10-0.15, we chose apparent magnitude limits for each cluster to select 500-600 candidate 



targets for a single Hectospec pointing (potentially up to 300 targets). The bright magnitude limit for these clusters 



is 17. 8< r <18.8, again corresponding to M* + 1 or fainter for most clusters (Figure 14) 



http: / / cas.sdss.org/ dr6/en/tools/chart /list. asp 
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Fig. 16. — Redshift versus projected radius for galaxies around HcCS clusters. The caustic pattern is evident as the trumpet-shaped 
regions with high density. The solid lines indicate our estimate of the location of the caustics in each cluster. Clusters are ordered 
left-to-right and top-to-bottom by decreasing X-ray luminosity. 
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Fig. 17. — Redshift versus projected radius for galaxies around HeCS clusters. The caustic pattern is evident as the trumpet-shaped 
regions with high density. The solid lines indicate our estimate of the location of the caustics in each cluster. Clusters are ordered 
left-to-right and top-to-bottom by decreasing X-ray luminosity. 
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Fig. 18. — Redshift versus projected radius for galaxies around HcCS clusters. The caustic pattern is evident as the trumpet-shaped 
regions with high density. The solid lines indicate our estimate of the location of the caustics in each cluster. Clusters are ordered 
left-to-right and top-to-bottom by decreasing X-ray luminosity. 
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Fig. 19. — Redshift versus projected radius for galaxies around HcCS clusters. The caustic pattern is evident as the trumpet-shaped 
regions with high density. The solid lines indicate our estimate of the location of the caustics in each cluster. Clusters are ordered 
left-to-right and top-to-bottom by decreasing X-ray luminosity. 
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Fig. 20. — Caustic mass profiles for the HeCS clusters. The thick solid lines show the caustic mass profiles and the thin lines show the lcr 
uncertainties in the mass profiles. The inner vertical solid line in each panel shows the radius T200- The next vertical line shows the smaller 
°f r 5.6 (the limit of bound structure) and r max (the maximum radius where the caustics are detected). For clusters with r max < rs.g, a 
dashed vertical line shows a lower limit on rs.g assuming no mass is present outside r max . Green shaded regions indicate the virial mass 
profile in the range (0.75-1.3)r200 (approximately from rsoo to rioo). 
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Fig. 21. — Caustic mass profiles for the HeCS clusters. The thick solid lines show the caustic mass profiles and the thin lines show the lcr 
uncertainties in the mass profiles. The inner vertical solid line in each panel shows the radius T200- The next vertical line shows the smaller 
°f r 5.6 (the limit of bound structure) and r max (the maximum radius where the caustics are detected). For clusters with r max < rs.g, a 
dashed vertical line shows a lower limit on rs.g assuming no mass is present outside r max . Green shaded regions indicate the virial mass 
profile in the range (0.75-1.3)r200 (approximately from rsoo to rioo). 
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Fig. 22. — Caustic mass profiles for the HeCS clusters. The thick solid lines show the caustic mass profiles and the thin lines show the lcr 
uncertainties in the mass profiles. The inner vertical solid line in each panel shows the radius T200- The next vertical line shows the smaller 
°f r 5.6 (the limit of bound structure) and r max (the maximum radius where the caustics are detected). For clusters with r max < rs.g, a 
dashed vertical line shows a lower limit on rs.g assuming no mass is present outside r max . Green shaded regions indicate the virial mass 
profile in the range (0.75-1.3)r200 (approximately from rsoo to rioo). 
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Fig. 23. — Caustic mass profiles for the HeCS clusters. The thick solid lines show the caustic mass profiles and the thin lines show the lcr 
uncertainties in the mass profiles. The inner vertical solid line in each panel shows the radius T200- The next vertical line shows the smaller 
°f r 5.6 (the limit of bound structure) and r max (the maximum radius where the caustics are detected). For clusters with r max < rs.g, a 
dashed vertical line shows a lower limit on rs.g assuming no mass is present outside r maI . Green shaded regions indicate the virial mass 
profile in the range (0.75-1.3)r200 (approximately from rsoo to rioo). 



